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Abstract 

The Jacobi polynomial has been advocated by many authors as a useful tool to 
evolve non-singlet structure functions to higher Q 2 . In this work, it is found that the 
convergence of the polynomial sum is not absolute, as there is always a small fluctuation 
present. Moreover, the convergence breaks down completely for large N . 

[PACS : 12.38.-t, 12.38.Bx, 12.90.+b, 13.90. +i ] 

The structure functions are the necessary tool in our effort to understand the hadronic 
structure and strong interaction. The study of Q 2 evolution of nucleon structure functions 
has been an important source of experimental information supporting Quantum Chromo- 
dynamics (QCD), which is believed to be the fundamental theory of strong interaction. It 
has already been shown some time back that [|T| QCD is the only theory which can explain 
the gross features of scaling violations in deep inelastic scattering (DIS). As a result a huge 
amount of effort is being put, both experimentally as well as theoretically, to understand the 
nucleon structure functions for different values of x and Q 2 . 

The evolution of quark distribution with Q 2 is goverened by the Altarelli-Parisi (AP) 
equation 0. To leading order in a s , the AP equation is given by, 



where q is the quark distribution, a s is the strong coupling, t = logQ 2 and P qq is the 

quark splitting function which represents the probability of a quark emitting a gluon and so 
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becoming a quark with momentum reduced by a fraction ~. Eq. ([!]), an integro-differential 
equation, is not very easy to solve; different methods have been proposed in the literature 
which can be grouped into three main categories. 

One method, already used by Altarelli, Nason and Rudolfi H is based on the assumption 
that, for a small variation of t, one can neglect the t dependence of the r.h.s in eq.(|l]) and 
realize, in steps, the evolution of quark distribution for a given St. 

The second method [Q consists of expanding the quark (parton) distribution functions 
into a truncated series of Chebyshev polynomial. This expanded form is substituted in eq. (|l|) 
and then the resulting coupled differential equations are solved self-consistently. 

The third method relies on the premise that the moments of the structure function (say 
F 2 ) depend only on Q 2 . This means that one can expand F 2 in terms of an ortho normal 
polynomial (usually Jacobi polynomial) such that the x dependence is carried by the polyno- 
mial whereas the full Q 2 dependence is confined to the weight factors. The variation of these 
weight factors with t can then be extracted from the knowledge of the variation of moments 
with t. In the present report we confine our attention to the use of Jacobi polynomial in 
solving the evolution equation of quark distrubution functions. 

The Jacobi polynomial was first used by Sourlas and Parisi || and later on elaborated by 
Barker et al. || and also Chyla et al. 0. We will follow the prescription of ref.0, though 
we have found that, both and M give similar results. For illustrative purposes, we restrict 
ourselves to the valence part of F 2 (x, Q 2 ). 

The method of orthogonal polynomials (here Jacobi Polynomial) is based on inverting 
moments with the help of orthogonal polynomials. The Jacobi polynomial is defined || as 

Qk a ' p (x) = J2 k CkAa,P)x J (2) 

3=0 

satisfying a weighted orthogonality relation, 

f'dxx^l - x) a Q k a ^(x)e l a ' /3 (x) = 8 kl (3) 
Jo 

Now the structure function F 2 can be expanded as 

F 2 (x, Q 2 ) = x%l - x) a Y,°°a k a %Q 2 )Q k a '%x) (4) 

k=0 
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From eq.([|) one can write the expansion coefficients a in terms of F 2 as 

a k al3 (Q 2 )= [ 1 F 2 (x,Q 2 )e k a ^(x)dx (5) 
Jo 

where the orthogonality relation of Qk a ^(x) (eq.^j has been made use of. 
Substituting eq.(Q) in the eq.([|) one gets, 

a k aP {Q 2 ) = Y, k C kj (*,P)vti + 2,Q 2 ) (6) 

3=0 

where the moments \i are given by, 

KJ,Q 2 )= f 1 dxx^ 2 F 2 (x,Q 2 ) (7) 
Jo 

We have used the general form of Jacobi polynomial |6|. |7J . 

Q k ^(x) = N k ^H k ^(x) (8) 
where the normalization factor is 

N k af3 = 9^(0) 
= (P+l 



(2k + a + (3 + l)r(Jfe + a + j3) 



T(k + a + l)T(k + /3+ l)k\ 
with (o) n = a(a + 1) • • • (a + n — 1) and ao = 1. The expansion coefficients C are then 



(9) 



The Jacobi Polynomial can be evaluated using eq.([|) or the recurrence relations for the 
polynomial^ @,|7|. 

In the following discussion we will study the reliability and convergence of Jacobi poly- 
nomial method for the evolution of structure functions. This is done firstly by starting from 
an analytical fitted formula for valence part of F 2 at a particular Q 2 and then reevaluating 
the F 2 at the same Q 2 using the Jacobi polynomials. We have used three different formulae 
at Q 2 = 3.5, 5 and 15 GeV 2 as given below 1,1,0, 

Q 2 = 3.5 GeV 2 

u(x, Qo 2 ) = y/x(l - x 2 ) 3 (0.594 + 0.461(1 - x 2 ) + 0.621(1 - x 2 ) 2 ) 

d( x , Qo 2 ) = y/x(l - x 2 ) 3 (0.072 + 0.206(1 - x 2 ) + 0.621(1 - x 2 ) 2 ) (11) 



g 2 = 5 GeV 2 



xu v (x, Q 2 ) = 1.78^(1 - x L51 ) 3 - 5 
xd v (x,Q 2 ) = 0.67x°' 4 (l - x L51 ) 4 - 5 

(12) 



Q 2 = 15 GeV 2 



xu v (x, Qq 2 ) — 777 n — pr^ a "(l — %) l3u 
B(a Ul p u + 1) 

xd v (x,Q 2 ) = — 1 —— x a *(l-xf« (13) 

-Dl"d, Pd + 1J 

where 5(a, (3) are the Euler beta functions, a u = 0.588±0.020±0.05, f3 u = 2.69±0.13±0.21, 
a d = 1.03 ± 0.10 ± 0.19 and /3 d = 6.87 ± 0.64 ± 0.80 ||. 

The variation of R = ^^(jormlia) w ^h the AT (number of terms of the Jacobi polynomial) 
for Q 2 = 15, 5 and 3.5 GeV 2 is shown in figure 1. For each Q 2 curves for x = 0.05, 0.4 and 
0.75 are plotted. We find that for Q 2 = 3.5 GeV 2 , there are large fluctuations in the R for 
smaller values of N. The value of R becomes 1 around N = 5 and stays the same upto a 
value of iV around 20. But immediately after that the value of R diverges. For larger values 
of x, R diverges for larger N. Similar features are present for higher values of Q 2 along 
with the additional feature that there are small oscillations in R around 1 till it diverges for 
N > 20. This behaviour of R shows that though there is an apprent convergence of F2 with 
N, this convergence may not be as conclusive as claimed by earlier authors 0, ■ 

The unreliability of the present method of evolution of structure functions is illustrated 
in figure 2. Here we have plotted the evolved value of only the valence part of from 
Qo 2 = 3.5 GeV 2 to Q 2 = 5 and 15 GeV 2 , for x = 0.05,0.4 and 0.75. Here we find that F 2 P 
diverges for iV > 12 which is much lower than the one observed (N > 20) in figure 1. Here 
again we find that the F2 P diverges earlier for lower values of x. The implication of these 
observations is that the even the apparent convergence of R does not say anything conclusive 
regarding the number of terms needed in the Jacobi Polynomial. 

In figures 3 and 4, we have plotted the evolved value of AF = F 2 P — F 2 n starting from 
Qq 2 = 3.5 and 5 GeV 2 respectively. The figures show that the AF gives same values for 
N < 14. For iV = 17 there is a large oscillation in the value of AF. The oscillations are 



larger for lower values of x and they decrease with increase in x. Furthermore, the evolved 
AF at Q 2 = 7, 12 and 15 GeV 2 have similar values. The above observations are valid for 
both Q 2 = 3.5 as well as Q 2 = 5 GeV 2 ; fluctuations being less for Q 2 = 5 GeV 2 . 

In conclusion we have studied the reliability and convergence of the method of Jacobi 
polynomials for the evolution of structure functions. We find that the convergence of this 
method is not wholly reliable in the sense that for large N, the valence part of F 2 diverges 
and there are large fluctuations in F 2 P — F 2 n . Furthermore, the final results at higher Q 2 
are strongly dependent on the initial fitted formula used, even for the value of N where 
fluctuations are not problematic. This is exemplified by the fact that even for N =10, the 
final value of F% — F£ at Q 2 = 15 GeV 2 is strongly dependent on whether one starts from 
Q 2 = 3.5 GeV 2 or 5 GeV 2 (figs. 3(c) and 4(c)). Thus, the method of Jacobi Polynomials, 
its simplicity and apparent success notwithstanding, is seen to be of limited validity in 
evaluating the Q 2 evolution of structure functions. It may be interesting and worthwhile to 
explore if there are any other polynomials which could perform better in this respect. 

We are grateful to X. Song for helpful discussions and to James S. McCarthy for helpful 
discussions as well as collaboration on a related project, the results of which would be 
published shortly. SKG would like to thank the Council of Scientific and Industrial Research 
(Govt, of India) for financial support. 
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FIGURE CAPTIONS 

Fig. 1. The ratio of the calculated value of the valence part of F$ and the corresponding 
value from the fitted formula (eq. (11-13)) , starting from a given value of Q 2 and reconstruct- 
ing it back at the same Q 2 using the Jacobi polynomial. Three curves in each plot correspond 
to x=0.05 (straight line), x=0.4 (long dash) and x=0.75 (short dash), respectively. 

Fig. 2. Evolved value of only the valence part of F$ (a) from 3.5 GeV 2 to 5 GeV 2 and 
(b) from 3.5 GeV 2 to 15 GeV 2 . 1, 2 and 3 in (a) and (b) correspond to x=0.05, 0.4 and 0.75 
respectively. 

Fig. 3. The evolved value of F 2 P - F 2 n from Q 2 = 3.5 GeV 2 to Q 2 = (a) 7 GeV 2 , (b) 12 
GeV 2 and (c) 15 GeV 2 . The evolution for N=6, 10 and 14 give the same value (solid line) 
whereas iV = 17 (dotted line) shows fluctuations for all the cases. 

Fig. 4. The evolved value of Ff - F 2 n from Q 2 = 5 GeV 2 to Q 2 = (a) 7 GeV 2 , (b) 12 
GeV 2 and (c) 15 GeV 2 . The evolution for N=Q, 10 and 14 give the same value (solid line) 
whereas N=17 (dotted line) shows fluctuations for all the cases. 
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